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ABSTRACT 

Colors of Trans Neptunian Objects (TNOs) are used to study the evolutionary processes of bodies in the out- 
skirts of the Solar System, and to test theories regarding their origin. Here I describe a search for serendipitous 
Sloan Digital Sky Survey (SDSS) observations of known TNOs and Centaurs. I present a catalog of SDSS 
photometry, colors and astrometry of 388 measurements of 42 outer Solar-System objects. I find a weak evi- 
dence, at the w 2(7 level (per trial), for a correlation between the g — r color and inclination of scattered disk 
objects and hot classical KBOs. I find a correlation between the g — r color and the angular momentum in the 
z direction of all the objects in this sample. These findings should be verified using larger samples of TNOs. 
Light curves as a function of phase angle are constructed for 13 objects. The steepness of the slopes of these 
light curves suggests that the coherent backscatter mechanism plays a major role in the reflectivity of outer 
Solar-System small objects at small phase angles. I find a weak evidence for an anti-correlation, significant 
at the 2(7 confidence level (per trial), between the g-band phase angle slope parameter and the semi-major 
axis, as well as the aphelion distance, of these objects (i.e., they show a more prominent "opposition effect" at 
smaller distances from the Sun). However, this plausible correlation should be verified using larger sample. I 
discuss the origin of this possible correlation and argue that if this correlation is real it probably indicates that 
"Sedna"-like objects have a different origin than other classes of TNOs. Finally, I identify several objects with 
large variability amplitudes. 

Subject headings: solar system; general — solar system; Kuiper Belt — techniques; photometric 



L INTRODUCTION 

Colors and variability of small bodies in the outer Solar 
System provide insight into the physical properties and evo- 
lution of these objects. The colors of these objects are be- 
lieved to be related to evolutionary processes such as colli- 
sions, resurfacing and the interaction of cosmic rays with the 
surface of these bodies (e.g.. Cooper et al. 2003; see how- 
ever Porter et al. 2010). The large body of color observa- 
tions of Trans Neptunian Object (TNOs; e.g., Luu & Jewitt 
1996; Delsanti et al. 2001; Hainaut & Delsanti 2002; Trujillo 
& Brown 2002; Tegler & Romanishin 2003; Almeida et al. 
2009; Santos-Sanz et al. 2009; Romanishin et al. 2010; Shep- 
pard 2010) is not entirely consistent with theoretical ideas (see 
the review in Jewitt, Morbidelli & Rauer 2008). The main 
characteristic of TNO colors is diversity. To date, the only se- 
cure correlation involving TNO colors is between the B — 
color and the orbital inclination of classical Kuiper Belt Ob- 
jects (KBOs) and scattered disk objectfl(e.g., Hainaut & Del- 
santi 2002; Trujillo & Brown 2002; Peixinho et al. 2008). I 
note that the reported correlations between the inclination and 
V — Rc or Rc — Ic colors remain controversial (e.g., Stephens 
et al. 2003). 

KBO variability is related to shape, binarity and albedo 
surface uniformity. Measuring the binary frequency allows 
of testing models of KBO binary formation (e.g. Goldreich, 
Lithwick & Sari 2002; Weidenschilling 2002). Moreover, in 
some cases binaries are used to determine masses (e.g., Noll 
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^ Here TNOs are defined as objects with semi-major axis larger than that 
of Neptune. KBOs and scattered disk objects are loosely defined and here we 
follow the definition of Morbidelli & Brown (2004). 



et al. 2004) and densities (e.g., Sheppard & Jewitt 2004; Gnat 
& Sari 2010) of KBOs. 

TNO variability studies typically require medium size tele- 
scopes and are therefore observationally demanding. The 
Sloan Digital Sky Survey (SDSS; York et al. 2000) provides 
imaging in the ugriz-hands of a considerable fraction of the 
celestial sphere. The photometric calibration of the SDSS is 
good to w 1% in the griz bands and w 2% in the u band (e.g.. 
Tucker et al. 2006; Padmanabhan et al. 2008). The SDSS as- 
trometric accuracy is w 0.1" (e.g.. Pier et al. 2003). However, 
given the short time interval within which the SDSS images 
were obtained (« 5 min), it does not allow in most cases the 
detection of KBO motiorQ Ivezic et al. (2001) and Juric et 
al. (2001) constructed a catalog of all the SDSS sources dis- 
playing a significant motion within the 5 min exposures — the 
SDSS Moving Object Catalog (SDSSMOC). However, in the 
fourth release of this catalog (SDSSMOC4) there are only 33 
entries of known objects with a > 10 AU. 

Here I describe a search for known small objects in the outer 
parts of the Solar System in the existing SDSS imaging data. 
A compilation of the photometric and astrometric properties 
of these bodies are presented and analyzed. The structure of 
this paper is as follows. In ij2]l describe the catalog of SDSS 
observations of outer Solar-System objects. In fj3]l discuss 
their colors, while in [2]I describe their variability properties. 
Finally, I summarize the results in fj5] 

2. A CATALOG OF SDSS OBSERVATIONS OF TNOS 

This section describes the construction of a catalog of SDSS 
observations of known outer Solar-System objects with semi- 
major axes a > 10 AU. 

^ At opposition, the typical geocentric on sky motion of a Solar-System 
object orbiting the Sun at 40 AU is ~ 0.25" in 5 min 

^ http://www.astro.washington.edu/users/ivezic/sdssmoc/sdssmoc.html 
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2.1. SDSS images whose footprints may contain known 

TNOs 

I retrieved a list of the orbital elements of all known (num- 
bered and unnumbered) minor planets in the Solar SysterrQ 
from the Jet Propulsion Laboratory (JPL) Horizon^ system 
(updated for 2010 August 2). Then, I selected all the objects 
with semi-major axis a > 10 AU. This list contains 1469 bod- 
ies. 

I used the SDSS (York et al. 2000) CasJobfl utiUty to gen- 
erate a catalog of all the imaged which are available in the 
SDSS database. Here, an "image" is defined uniquely by the 
SDSS run, rerun, camcol (camera column), and fielcF^ The 
catalog contains all the images included in the SDSS data re- 
lease 7 (DR7), Segue, and Stripe 82 (Abazajian et al. 2009). 
For each image 1 obtained the time at which it was observed 
and I calculated the coordinates of its four cornerfl. 

Next, 1 used the JPL Horizons system to generate daily 
ephemerides for each of the 1469 objects between Julian 
Day (JD) 2451070 and 2454467. This JD range contains 
all the SDSS observations in DR7. For each entry in the 
daily ephemerides of each object 1 checked whether it is con- 
tained within any of the polygons defining the corners of all 
the SDSS images taken within one day of the ephemeris en- 
try. If a match was found, then the object ephemeris was re- 
generated for the exact time at which the image was taken (to 
an accuracy of 1 min). In total, 4642 possible observations of 
574 outer Solar-System objects were found. Of these 845 en- 
tries are of objects with a predicted V-band magnitude, at the 
time of observation, brighter than Vpred — 22 mag. 

2.2. Photometry and astrometry of TNOs in SDSS images 

Next, I searched for sources in the SDSS images near the 
predicted position of the outer Solar-System objects. Unlike 
"typical" minor planet surveys, this method only yields a sin- 
gle image per object per field, so one cannot use the motion 
of the object between two images of the same field to ver- 
ify whether it is indeed a Solar-System object (rather then a 
variable star or a transient). Therefore, as described below, 1 
exercised great care to remove false identifications or contam- 
ination by nearby sources. 

For each entry in the catalog of SDSS images possibly con- 
taining an observation of a Solar-System obje ct w ith a pre- 
dicted magnitude brighter than V^red ~ 22.0 ( i]2.1l ). 1 down- 
loaded the SDSS source catalog corresponding to that im- 
agj^. Then 1 searched for all the SDSS sources within 8" 
of the predicted position of the object. Figure [T] shows the 
distribution of the angular separations between the predicted 
object position and the nearest SDSS source. 

In order to remove photometric measurements contami- 
nated by nearby sources and possible false detections, 1 se- 
lected only sources which have exactly one SDSS match 
within 8" of the outer Solar-System object predicted position. 
1 further demanded that this SDSS source is within 1 .5" from 
the predicted position of the object. Moreover, I selected only 

^ http://ssd.jpl. nasa.gov/?sb_lem 
^ http://ssd.jpl.nasa.gov/7horizons 
' http://casjobs.sdss.org/casjobs/ 

1,589,702 images. 

" see definitions in: http://www.sdss.org/dr7/glossary/index.html] 

Performed by transforming the SDSS great circle coordinates to equa- 
torial coordinates; http://www.sdss.org/dr7/products/general/astrometry.html 
" These are the tsObj files stored in http://das.sdss.org/imaging/ which are 
described in http://www.sdss.Org/dr7.l/dm/flatHles/tsObj.html 
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Fig. 1 . — Histogram of angular distances between the predicted position of 
the Solar-System objects and the position of the SDSS source nearest to the 
predicted position (either the Solar-System object or a background object). 
This plot is shown for 672 measurements in which I find at least one SDSS 
source within 8" of the predicted position of the object. 



sources for which there are no USNO-Bl (Monet et al. 2003) 
objects within 8" from the position of the SDSS matched 
source. Finally, 1 removed from the list of SDSS matches 
entries with r-band magnitude errors larger than 0.2 mag. The 
final list contains 388 probable identifications of 42 unique 
outer Solar-System objects in SDSS images. Each of these 
objects has between 1 and 49 measurements. Table [Upresents 
the astrometric properties of the 388 identifications, while Ta- 
ble |2] lists all the predicted and measured photometry. 

I also calculated the absolute planetary magnitud43. ne- 
glecting phase effects (see SjH, Hf — m/ — 51ogio(7?A), where 
m/ is the apparent magnitude in band / (m, g, r, i or z), R 

is the object's heliocentric distanc43, and A is its geocentric 
distance. The values of R, A and the phase angle j3 (defined 
as the Sun-target-observer angle), for each observation, are 
listed in Table |2l 1 also calculated the median, standard de- 
viation (StD), and range of the absolute magnitude measure- 
ments (Table 12]). 

2.3. Verification 

As shown in Figure[Tl a large fraction of the matched SDSS 
sources are found within 1 .5" from the predicted position of 
the Solar-System objects. In contrast, the probability for a 
match with background sources should increase as the square 
of the distance. Since the number of matches above 1.5" is 
small 1 argue that the fraction of false identification and con- 
taminated photometry in Tables [T}{2] is negligible. Further- 
more, 1 note that assuming a source density of lO'^^deg^^ in 
the SDSS images, the probability of finding a source within 
1 .5" from a random position is 0.5%. 

Nevertheless, as an additional test 1 uploaded cutouts of 
the SDSS images containing some of the candidates, along 
with SDSS images of the same sky positions taken at differ- 
ent epochs. If such an extra epoch image was not available 
then 1 uploaded instead an image from the Palomar Sky Sur- 

''^ Defined as the magnitude of an object observed at opposition and at 
1 AU from the Sun and Earth. 

Denoted by R to distinguish it from the SDSS r-band magnitude. 



KBOs observations in SDSS 
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TABLE 1 

ASTROMETRIC MEASUREMENTS OF SOLAR-SYSTEM OBJECTS IDENTIFIED IN SDSS IMAGES 



Name 


JD-245()000 
day 


Run 


Rerun 


Camcol 


Field 


"2000 
deg 


sobs 
°2000 
deg 


"2000 
deg 


c-pred 
"2000 

deg 


Dist 
arcsec 


24835 


1464.856934 


1035 


40 


4 


126 


24.909451 




24.909542 


14.267833 


0.43 


26375 


2318.829706 


2986 


40 


3 


273 


159.458594 


6.417053 


159.458667 


6.417111 


0.34 




2338.709340 


3015 


40 


3 


300 


159.008464 


6.578481 


159.008458 


6.578528 


0.17 


35671 


1819.838123 


1755 


40 


2 


340 


354.611093 


-0.825181 


354.611083 


-0.825139 


0.16 




2196.787392 


2649 


40 


2 


94 


355.871559 


-0.418224 


355.871583 


-0.418222 


0.09 



Note. — Astrometric measurements of 388 identifications of 42 outer Solar-System objects. Name is the minor planet number or provisional designation, 
Run/Rerun/Camcol/Field identifies the unique SDSS ugriz image, while JD provides the time at which the r-band image was obtained, a and 5 are the J2000.0 
coordinates of the object, while superscript "obs" is for observed coordinates, and superscript "pred' is for predicted coordinates. Dist is the angular distance 
between the observed and predicted coordinates. This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion of the full 
table is shown here for guidance regarding its form and content. 



TABLE 2 

Photometric properties of Solar-System objects identified in SDSS images 



Name 


JD-2450000 


R 


A 


/3 


Vprcd 


Type u 


S 


r 


i 


z 


A« 




Ar 


Ai 


Az 




day 


AU 


AU 


deg 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


24835 


1464.856934 


39.497 


38.512 


0.24 


20.8 


6 


20.879 


20.363 


20.192 


19.930 


0.271 


0.038 


0.034 


0.043 


0.145 


26375 


2318.829706 


34.179 


33.221 


0.42 


20.1 


6 


20.859 


20.087 


19.755 


19.634 


0.269 


0.033 


0.025 


0.026 


0.059 




2338.709340 


34.191 


33.206 


0.18 


20.0 


6 


20.929 


20.161 


19.784 


19.592 


0.605 


0.034 


0.027 


0.026 


0.064 


35671 


1819.838123 


38.156 


37.188 


0.38 


21.6 


6 


21.760 


21.262 


21.140 




0.425 


0.076 


0.062 


0.074 


0.275 




2196.787392 


38.109 


37.203 


0.63 


21.7 


6 


21.653 


21.270 


21.291 




1.087 


0.063 


0.057 


0.079 


0.242 



Note. — As in Table [T] except for the photometric properties. Vpi-ej is the minor planet predicted V-band magnitude at the time of the SDSS observation. 
Type is the SDSS morphological classification (6: star; 3: galaxy), ugriz are the SDSS magnitudes, while their corresponding errors ai'e Ah, Ag, Ar, Ai and Az. In 
cases in which the error magnitudes are larger than 0.2 mag, 1 replaced the SDSS magnitude with the "no-data" symbol. However, I kept the errors in the table. 
Note that the absolute planetary r-band magnitude of measurement number 18 of object 145452, measurement number 7 of object 145453, and measurement 
number 13 of object 145480 deviate by more than one magnitude from the median absolute planetary magnitude and are probably bad measurements. These 
measurements are listed in this table but are not used in the subsequent analysis (e.g., they are not shown in Figure[3]and they are excluded from the phase-angle 
slope parameter fits summarized in Table O. This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion of the full 
table is shown here for guidance regarding its form and content. 

ve\P1 (Reid et al. 1991). I inspected by eye about 100 of these 
cutouts and verified that the Solar-System object candidate is 
indeed detected only in one epoch. Finally, 1 note that the 
differences between the predicted and measured magnitudes 
(Table|2]l istypically small, on the order of 0.3 mag. 

In TableUl list also the SDSS morphological type (6: star; 
3: galaxy) for the TNOs. Some of the sources are identified 
as possible resolved objects. This is presumably because the 
reliability of SDSS star-galaxy separation degrades near the 
survey detection limit. 

3. OBJECT COLOR 

The absolute planetary magnitude of objects identified in 
SDSS images are listed in Table |3] If multiple-epoch ob- 
servations are available, 1 adopt the median of the object 
magnitudes over all epochs as the object's magnitude. In 
cases in which an object was observed in multiple epochs, 
I also give the standard deviation (StD) of the absolute mag- 
nitudes, and the range of the r-band absolute magnitudes. 1 
note that the variability indicators in this table do not sep- 
arate between variability due to phase-angle variations and 
other causes (e.g., rotation). Separation of phase-angle and 
rotation-induced variability is possible only when a large 
number of observations is available. Nevertheless, the vari- 
ability indicators in this table give a rough idea regarding 
which objects may be variable and which objects are less 
likely to be variable. For objects which have more than ten ob- 
servations more reliable variability indicators, which are cal- 

http://archive.stsci.edu/cgi-bin/dss_form 



culated after subtracting the phase-angle variations, are avail- 
able in Table 5. Figure |2] shows the § — r vs. r — i color-color 
diagram for the 37 objects for which both g — r and r — i color 
measurements are available. The symbols indicate different 
subclasses of objects (see figure legend and caption). I note 
that the § — r vs. r — i locus of objects in Figure|2]is generally 
similar, but offset, relative to the B — V vs. R I colors locus 
of TNOs (for comparison see Fig. 2 in Tegler & Romanishin 
2003). 

In order to explore possible correlations of the colors with 
orbital properties I calculated the Spearman rank correlation 
coefficients between the g — r, r — i, and g — i colors and var- 
ious photometric and orbital properties of these objects and 
their subgroups (e.g., the groups listed in Fig. |2]i. This ap- 
proach has the disadvantage that it introduce many trials, and 
reduces the significance of any reported correlation by a com- 
plicated way that depends on the correlations between the var- 
ious parameters. Nevertheless, this may give us some ideas 
about where to look for correlations when larger samples, 
based on the same filters, become available. For each correla- 
tion 1 also calculated the probability of getting a value larger 
than the correlation coefficient. This was calculated from the 
correlation coefficients' probability distributions derived from 
10^ bootstrap simulations (Efron 1982; Efron & Tibshirani 
1993). In each simulation, considering two lists "X" and "F", 
1 select for each entry in X a random entry in Y and calculate 
the correlation between the two randomly permuted lists. I 
calculated the correlations of the three colors g — r, r — i, and 
g — / with the r — / color, Hr, semi-major axis a, orbital eccen- 
tricity e, orbital inclination /, perihelion distance q, aphelion 
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TABLE 3 

Mean photometric properties of outer Solar-System objects identified in SDSS images 



Name 




Hu 
mag 


mag 


Hr 

mag 


H, 
mag 


Hz 
mag 


USlD 

mag 


gSiD 

mag 


rsiD 
mag 


iSTD 

mag 


ZStD 

mag 


^ range 

mag 


a 
AU 




e 


/ 

deg 


24835 


1 




4.968 


4.452 


1 

4.281 


4.019 




fj 
u 




A 
U 


nnn 


n 


nnn 

UUU 


n 
u 


nnn 

UUU 


n 

y) 


nnn 

UUU 


A 1 


QS7 





106 


97 
Z / 


nnn 

l/UU 


26375 


2 




5.619 


4.848 


4.494 


4.337 




Q 


050 


Q 


052 


Q 


021 


Q 


029 


Q 


073 


J J 


108 





415 


7 


D-? i 


35671 


16 




6.079 


5.573 


5.450 


4.975 






104 


Q 


106 


() 


nQQ 

yjyy 


u 


028 


() 


398 


JO 


1 10 





043 


4 


596 


3»62s (Huya) 


1 




5.600 


4.775 


4.394 


4.249 




u 


nnn 


n 
u 


nnn 


n 


nnn 

UUU 


n 
u 


nnn 

UUU 


A 
U 


nnn 

UUU 


jy 


3 13 





276 


1 s 


ASS 


D54»y (Leto) 


1 




7.182 


6.519 


6.136 


5.893 




fJ 
u 


nnn 
uuu 


A 

u 


nnn 


n 
u 


nnn 

UUU 


n 
u 


nnn 

UUU 


n 

y) 


nnn 

UUU 


QQ 

yy 


D /O 





821 


99 
ZZ 


-jZJ 


73480 


1 1 




9.270 


8.641 


8.337 


8. 15 1 






1 OS 


A 

yj 


U /*+ 


A 

u 


n7Q 

y) ly 


n 
u 




n 
u 


ZJO 


J 1 


9ilS 
Z'tJ 





572 


1 fi 
ID 


f^97 
DZ / 


79360 


4 




6.187 


5.287 


4.773 


■ ■ ■ 




A 

u 


1 'iO 


A 
U 


1 


n 
u 


Uj'-f 






n 
u 


^9^ 

jZj 




1 jy 





008 


z 


9'^n 

zjU 


82075 


3 




5.222 


4.496 


4.215 


4.084 




n 
u 


1 S'^ 


A 
Lf 


1 1 Q 


A 

u 


UDU 


n 
u 


1 QQ 

1 yy 


A 

u 


99S 
ZZJ 


S7 


DoD 





288 


1 Q 

ly 


s4n 

o'+U 


82155 


1 




6.634 


6.141 


5.693 








nnn 


A 
U 


nnn 

UUU 


n 
u 


nnn 

UUU 






n 
u 


nnn 

UUU 


9.A 


OZJ 





617 


1 9 
iZ 


710 

/ 3y 


82158 


1 




6.963 


6.129 


5.708 






A 

u 


uuu 


A 
U 


nnn 
uuu 


n 
u 


nnn 

UUU 






n 
u 


nnn 

UUU 


919 
Z iZ 


fSOO 





839 


in 

jU 


9.A1 
o4Z 


90482 (Orcus) 


2 


4.110 


2.702 


2.222 


2.082 


2.054 


0.021 


A 

u 


n^n 


A 
U 


n7n 
u / u 


u 


ni fi 

Ui o 


A 

u 


n7^ 

U / J 


u 


nQQ 

Kiyo 


TQ 

jy 


1 71 





227 


9n 
zU 


^71 

3 13 


1 19878 


2 




7.120 


6.216 


5.674 


5.198 




A 

\) 


Id J 


n 


U'-FU 


A 


UM-O 


n 
u 


nnn 

UUU 


n 


UJD 


^'K 

J J 


D^-O 





344 


1 s 


7SQ 

1 3y 


120132 


1 




5.436 


4.653 


4.265 


4.331 







000 


Q 


000 





000 


Q 


000 


Q 


000 


49 


1 Q7 

ly 1 





247 


1 1 


7QS 

/ yo 


120181 


2 




8.118 


7.309 


6.845 








043 


Q 


075 


() 


043 









106 


32 


J JU 





111 


2 


111 


135182 


1 


















Q 


000 













000 


37 


285 





021 


]^ 


837 


1 mnn z 
ijy/ ID 


10 




7.806 


7.047 


O.JV3 






A 
yj 


UD / 


n 
u 


1 Q7 


A 
U 








n 




'?Q 

jy 


yiHy 





199 


o 


ASn 

'foU 


1/1/1 QC\n 


17 




5. 148 


4.315 


4.026 


3.830 




A 

y) 


UoO 


n 


UDV 


n 


Z.DZ. 


n 
u 


1 79 
1 /z 


{\ 

y) 


99"^ 
LI. J 


^Q 

jy 


9ni 

ZUj 





043 


Q 

y 


JZt 


145451 


27 


5.996 


4.852 


4.436 


4.289 


4.283 


0.251 


A 

\) 


nsA 


A 

u 


n7n 
u /u 


n 
u 


nsA 


n 
u 


1 9A 
1 zo 


n 
u 


^9 1 

jZl 


Q 1 

yV 


79n 
/ zu 





617 


19. 
Zo 


7SQ 

1 Dy 


145452 


49 


5.922 


4.489 


3.694 


3.370 


3.206 


0.000 


A 
V) 


nQ9 


A 

u 


n7S 

U / J 


u 


UD / 


n 
u 


non 

yjyyj 


A 
U 


^7S 


41 


7SQ 

/ Jy 





028 


1 Q 

ly 


91f^ 
ZjD 


145453 


36 


5.740 


4.454 


4.003 


3.898 


3.843 


0.331 


A 
yj 


UD / 


A 
U 


n'^7 

U J / 


u 


UJ-J 


n 
u 


1 AS 


A 
U 


Ml 


^-j 


All 

"rLZ. 





143 


19 
Zo 


sno 

JU7 


1 A^AQCX 

i4j4oU 


40 




5.235 


4.45 1 


4. 1 19 


3.722 









Q 




u 


100 


Q 


1 SA 


Q 


852 


1 yj 


Jyv 





397 


26 


19 Q 


150642 


2 




6.463 


5.935 


5.738 


■ 




Q 




Q 




Q 


409 









202 


H-J 


ni Q 
uiy 


A 

u 


1 Id 


10 


Z J J 








4.221 


3.451 


3.093 


2.936 




A 




A 
U 


nnn 

UUU 


n 

y) 


nnn 

UUU 


n 
u 


nnn 

UUU 


n 
u 


nnn 

UUU 
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Note. — Columns description: N^^i^^ is the number of r-band observations in Table 121 H„ through are the median absolute planetary magnitudes, not 
coiTected for phase angle, in the ugriz bands; u,,^ through Zstd are the StD in //„ through H-, respectively, but after removing the three bad measurements (see 
Table|2j; and r^ange is the range in Hy over all measurements, excluding the three bad measurements. 



distance Q, the orbital angular momentum in the z-direction 
(L^ = yja{\ — e^^cos/), and the Tisserand parameter calcu- 
lated with respect to Neptune {Jn = '-^ +2^J^^{l~e^cosI, 

where = 30. 104 AU is the orbital semi-major axis of Nep- 
tune). 

The only significant correlation reported in the literature is 
the color-inclination correlation (e.g., Trujillo & Brown; Peix- 
inheo et al. 2008). Here, I find only weak evidence for this 
correlation. Specifically, I find that the g — r colors have a 
correlation coefficient of —0.90 and —0.68 for scattered disk 
objects and hot KBOs, respectively. The probability of get- 
ting correlation coefficient which are smaller than these val- 
ues are 1.4% and 2.3% (per trial), respectively. For the rest 
of the populations investigated here the g — r-inclination cor- 
relation one-sided false alarm probability is larger than 2.5% 
(corresponds to 2a). The correlation I find is weaker and less 
significant than that found in other studies. Possible explana- 
tions for the differences between the correlations found in this 
paper and in other works are: (i) the different filters used by 



different studies; (ii) selection biases that plague the various 
samples; and (iii) the small sample size. 

The nature of the color-inclination correlation is not clear 
Among the possible explanations are collisional resurfacing 
(Luu & Jewitt 1996; Jewitt & Luu 2001) in which collisions 
between TNOs expose fresh material and change their colors 
and at the same time excite their inclinations. Another possi- 
bility is that the colors of KBOs are primordial and related to 
dynamical groupings. However, both explanation have been 
challenged by observations (see Trujillo & Brown 2002; Volk 
&Malhotra 2011). 

Peixinho et al. (2008) found that there is a break in the "re- 
lation" between color and orbital inclination, where objects 
with / 12 deg shows no correlation with color. Moreover, 
the perihelion distance (= fl[l — e]) and inclination of Classi- 
cal KBOs are known to be loosely related. I note that the cor- 
relation of color with inclination and perihelion distance have 
a functional resemblance to the functional form of the orbital 
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Fig. 2. — The g — r vs. r — i color-color diagram for 37 objects. Different 
symbols represent distinct classes of objects. Here, I define resonant objects 
as those having orbital periods within 2.5% of 3/2, 4/3, 2, or 5/2 times the 
Neptune orbital period. Scattered disk objects are defined as those having 
a > 28 AU, q < 32 AU, g > 32AU and that are not resonant objects. Cold 
classical objects are defined as those having a > 30 AU, a < 60 AU, / < 6 deg 
and not being resonant or scattered disk objects. Hot classical are defined 
similarly to cold classical but having / > 6 deg. Centaurus objects are defined 
as those having 28 AU > a > 5.204 AU, and "Sedna'-like objects are all the 
objects with a > 60 AU that are neither resonant objects nor scattered disk 
objects. I note that some of the various subclasses are not well defined. 

angular momentum in the z-direction Curiously, I find 

that the correlation between the g — r color and L, has a one- 
sided false alarm probability of 0.4%. It will be interesting to 
test this correlation using larger samples. 

I note that collisions between objects conserve the total an- 
gular momentum of the bodies involve in the collision. Since 
the angular momentum of individual bodies is not conserved, 
the absolute value of the angular momentum of individual 
bodies may (at least for mostly elastic collisions) statistically 
increase after a collision. However, this depends on the details 
of the collisions (e.g., elasticity and initial orbits). Therefore, 
I cannot rule out that the Lz-color correlation, if real, is a 
byproduct of TNO collisions. At this stage, it is not clear that 
the color-L- correlation really has a physical meaning rather 
than being a combination of several (physically unrelated) 
correlations (e.g., color-inclination, and perihelion distance- 
inclination correlations). 

4. VARIABILITY 

Thirteen objects in my sample have ten or more SDSS r- 
band measurements with errors smaller than 0.2 mag. Al- 
though the observations are too sparse to unambiguously 
identify periods, they are good enough to study the objects' 
reflectivity as a function of phase angle ( H4.li and to search 
f or la rge amplitude variability due to rotation and binarity 

(E3. 

4. 1 . Phase angle variations 

Solar-System bodies are known to vary in brightness with 
phase angle. There are two important physical reasons for 
this variation. The first is shadow hiding, in which particles 
on the planetary surface cast shadows on adjacent areas: the 
shaded area is minimized near opposition. The second is an 

L- depends on cos(/) which varies by only 2% between and 12 deg. 



interference mechanism called coherent backscatter in which 
reflected light, depending on the regolith properties, may con- 
structively interfere, resulting in an increased brightness at op- 
position (Hapke 1993; 2002). 

Hapke (2002) presented models of these effects. These 
models have seven degrees of freedom. Given the relatively 
small number of observations and limited range of phase an- 
gles in which the SDSS observations were obtained, I fit a 
linear relation of the form 



Hfip)=Hf^o + Sfli. 



(1) 



Here, / is the filter name (g, r or /), Hf Q is the absolute plan- 
etary magnitude at zero phase angle, 5/ is a linearized phase 
angle slope parameter for filter /, and j3 is the phase angle. 
Figure [3] shows as a function of j3 for these 13 objects. 

As seen in other Solar-System objects, these 13 objects with 
one exception, are brightest near opposition. To quantify this, 
in Table |4] I summarize the phase angle slope parameters and 
related information for the 13 objects. The only object that 
does not follow this rule is 139775. A plausible explanation 
is that this object has a la rge intrinsic variability due to rota- 
tion or binarity (see i l4.2b . I note that the fits are performed 
only for measurements which //r(j3) is within 1 magnitude of 
the median of //, (j3). The three measurements which do not 
fulfill this condition, and were removed, are indicated in Ta- 
ble |2] Interestingly, in most cases the slope parameter 
) is larger than 0.04 mag deg^ ' . Such large slope param- 
eters were argued to be the result of coherent backscatter (see 
Schaefer, Rabinowitz & TourteUotte 2009). 

Although the sample of objects for which I measure the 
slope parameters is small I attempted to look for correlations 
between the slope parameters and the orbital parameters of 
these TNOs. Any hints for correlations found here can be 
tested in the future using larger samples. The only notable 
anti-correlations I find are between Sg and a, Q, and P, where 
P is the orbital period. For example, the correlation between 
and a is —0.59, and the probability to get a correlation 
smaller than this is 1.8% per trial (roughly 2(7 significance). 
In order to test if this correlation is real, larger samples are 
required. Figure |4] presents Sg as a function of a. This figure 
suggests that most of the apparent correlation arise due to a 
difference between objects with large aphelion distances (i.e., 
"Sedna"-like orbits) and the rest of the population. 

This anti-correlation means that the g-band slope parameter 
is steeper for objects which are closer to the Sun. However, 
this finding is based on a small sample of only 13 objects. 
A possible selection bias that may cause such a correlation 
is that closer objects are visible also on larger phase angles, 
whereas further objects are visible only near j3 « 0. Since 
in reality the derivative of the absolute-magnitude phase- 
angle relation increases (in absolute value) near opposition, 
this may introduce the observed anti-correlation between the 
slope parameter and the semi-major axis. However, most 
of the SDSS observations were taken near a phase angle of 
j3 = 1 deg (see Figure O. I also note that this selection bias 
should mostly induce a correlation with q rather than with a 
or Q, since objects with smaller q are easier to detect. Another 
possible caveat is that for three objects I find negative g-band 
slope parameters (see Table |4|i presumably affected by mea- 
surement errors and/or variability. Therefore, more observa- 
tions are required in order to confirm the existence of such a 
correlation. 

If this correlation is real, then there are several possible ex- 
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Fig. 3. — Absolute planetary r-band magnitude as a function of phase angle (/3). The negative slope parameter of 139775 is presumably due to variability. 



TABLE 4 

Linearized phase angle slope parameters for objects with multiple observations 



Name 


Hr.O 

mag 


5,- 

mag deg^ ' 




dofr 


mag 


mag deg ' 




dofg 


35671 


5.43 ±0.04 


0.215±0.053 


10.6 


14 


6.02 ±0.04 


0.111±0.058 


27.5 


14 


73480 


8.21±0.15 


0.122 ±0.043 


21.6 


9 


8.83±0.16 


0.128 ±0.046 


28.4 


8 


139775 


7.08 ±0.07 


-0.161 ±0.092 


20.0 


8 


7.66±0.30 


0.224 ±0.484 


0.0 





144897 


4.25 ±0.01 


0.141 ±0.023 


28.4 


15 


5.03 ±0.02 


0.195 ±0.030 


29.4 


15 


145451 


4.31 ±0.02 


0.171 ±0.020 


61.9 


25 


4.78 ±0.02 


0.096 ±0.022 


121.6 


25 


145452 


3.50±0.01 


0.160 ±0.011 


91.1 


46 


4.26 ±0.02 


0.184±0.017 


108.2 


45 


145453 


3.98±0.01 


0.031 ±0.015 


65.6 


33 


4.49 ±0.01 


-0.038±0.017 


125.1 


33 


145480 


4.26 ±0.08 


0.157 ±0.078 


70.2 


37 


5.24±0.12 


-0.035 ±0.107 


59.8 


32 


2002 QX47 


8.69 ±0.05 


0.081 ±0.029 


15.1 


10 


9.20 ±0.06 


0.116±0.036 


9.1 


9 


2003 QW90 


4.96 ±0.06 


0.098 ±0.085 


14.4 


9 


5.92 ±0.09 


-0.028 ±0.142 


5.9 


6 


2004PG115 


4.94 ±0.07 


0.081 ±0.054 


48.4 


11 


6.05 ±0.13 


-0.1 12 ±0.096 


7.0 


10 


2005 R043 


7.03 ±0.05 


0.222 ±0.056 


18.9 


30 


7.66 ±0.06 


0.155 ±0.067 


40.5 


24 


2005 RS43 


4.94 ±0.02 


0.163 ±0.038 


57.6 


45 


5.60 ±0.03 


0.126 ±0.047 


76.3 


38 


Note. — 


Mean photometric properties and slope parameters for the 13 outer Solar-System objects with 


more than nine SDSS observations. 


The fits are 



performed only using measurements with photometric errors smaller than 0.2 mag. Therefore, the g-band and r-band slope measurements are not always based 
on data points taken at the same epochs. Column descriptions: Hf Q is the best fit absolute planetary magnitude at zero phase angle for filter /; Sf is the linearized 
slope parameter; xj and doff indicate the and the number of degrees of freedom of the best fit for filter /. 
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Fig. 4. — The linearized g-band slope parameter 5^ as a function of the 
semi-major axis a of the 13 objects with more than nine photometric mea- 
surements. Symbols as in Fig. |2] The typical errors in Sg are roughly 
0.05 mag deg"' (see Table|4). 

planations: (i) the regolith or surface composition properties 
of TNOs vary with distance from the Sun; or (ii) "Sedna"-Hke 
objects have distinctive surface properties which are related 
to their origin. The variation in surface properties as a func- 
tion of distance from the Sun can originate, for example, if 
there are variations in the impact rate with micro-meteoroids 
as a function of heliocentric distance, or due to the crystal- 
lization properties of some ices. Based on the Voyager-I and 
II spacecraft measurements, Gurnett et al. (2005) argued that 
the number density of dust particles, as a function of helio- 
centric distance, is roughly uniform (up to distance of about 
100 AU). Since the typical orbital speed of objects at 90 AU 
is 0.6 of that of objects at 30 AU from the Sun, this implies 
that the micro-meteoroids impact rate does not change dra- 
matically for objects in my sample. 

As different ices freeze at different temperatures, the sur- 
face properties may also be affected by the equilibrium tem- 
peratur43 and escape velocity from the object (e.g., Schaller 
& Brown 2007). For an albedo of A = 0.04, the equilibrium 
temperature varies between about 50 K at 30 AU from the Sun 
to 29 K at a heliocentric distance of about 90 AU. However, 
the three "Sedna"-like objects shown in the right-hand side of 
Figure Hare 145451, 145480 and 2004 PG115. During the 
SDSS observations, these objects were near perihelion at dis- 
tances of 35, 46 and 36 AU from the Sun, respectively. There- 
fore, their actual equilibrium temperature, at the time of ob- 
servations, were similar to those of some of the other objects 
in Figured 

I conclude that differences in surface properties induced by 
the current orbit of these objects are unlikely. However, I can- 
not rule out that "Sedna"-like objects have a different origin 
than some of the other classes of TNOs. 

I note that Schaefer et al. (2009) found a significant cor- 
relation between the slope parameter (near phase angle of 
j3 = 1 deg) and the B~ color index and also a possible ex- 
cluded region in the slope parameter vs. inchnation phase 



space. I do not find any indication for such correlations. 
However, the Schaefer et al. (2009) sample is larger (35 ob- 
jects) and contains more diverse planetary objects including 
the largest KBOs. 

4.2. Variability due to rotation and binarity 

Tableland Figure|3]suggest that some of the objects in the 
catalog presented here are variable. The variability of these 
sources may be the result of one or more of the following rea- 
sons: Small minor planets, probably with radii smaller than 
~ 100 km, may have irregular shapes and therefore are vari- 
able. Moreover, TNOs may show variations due to inhomo- 
geneous surfaces. Alternatively, fast rotation of objects held 
by their own gravity (i.e., radii larger than ^ 100km) may 
induce a highly non-spherical equilibrium configuration and 
therefore large amplitude variations (e.g., Leone et al. 1984; 
Rabinowitz et al. 2006). Finally, contact binaries may show 
prominent eclipses (e.g., Sheppard & Jewitt 2004; Gnat & 
Sari 2010). 

Objects of 100 km radius with an albedo of 0.04 (0.1) will 
have absolute planetary magnitude of 7.6 (6.6). Therefore, 
most of the objects in Table|4|are probably larger than 100 km. 
In this case, it is probable that large amplitude variations are 
either due to fast rotation, inhomogeneous surface albedo, 
or binarity. Since all these possibilities are interesting, pho- 
tometric follow-up observations of the most highly variable 
sources in Tables [3] and |4] are desirable. 

5. SUMMARY 

I cross-correlate SDSS observations with the ephemerides 
of Solar-System bodies with a > 10 AU. I present a catalog 
of SDSS photometric and astrometric measurements of such 
minor planets based on SDSS observations. After removing 
possible contaminated measurements. I am left with 388 ob- 
servations of 42 unique objects. 

I find weak evidence for the previously reported inclination- 
color correlation in the scattered disk objects and hot classi- 
cal KBOs. I find marginally stronger correlation between the 
g — r color and orbital angular momentum in the z direction, 
L^, of the entire population studied here. I note that a corre- 
lation with L, is consistent with the finding of Peixinho et al. 
(2008) that objects with inclination below about 12 deg shows 
no color-inclination correlation. 

Finally, the method presented here to collect photometric 
observations of minor planets in surveys which were not de- 
signed for Solar-System observations can be utilized in other 
ongoing and planned surveys such as the Palomar Transient 
Factory (Law et al. 2009; Rau et al. 2009), Pan-STARRS 
(Kaiser et al. 2002), SkyMapper (Keller et al. 2007), and 
LSST (Tyson et al. 2003). 
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